arXiv: 1508.04461 vl [physics.class-ph] 17 Jul 2015 


Wave propagation in complex coordinates 

S. A. R. Horsley/ C. G. King/ and T. G. Philbin^ 

^Department of Physics and Astronomy, 

University of Exeter, Stocker Road, Exeter, EXf 4Q^ 

Abstract 

We investigate the analytic continuation of wave equations into the complex position plane. For 
the particular case of electromagnetic waves we provide a physical meaning for such an analytic 
continuation in terms of a family of closely related inhomogeneous media. For bounded permittivity 
profiles we find the phenomenon of reflection can be related to branch cuts in the wave that 
originate from poles of e{z) at complex positions. Demanding that these branch cuts disappear, we 
derive a large family of inhomogeneous media that are reflectionless for a single angle of incidence. 
Extending this property to all angles of incidence leads us to a generalized form of the Poschl Teller 
potentials. We conclude by analyzing our findings within the phase integral (WKB) method. 
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The propagation of linear waves through inhomogeneous media is straightforward to com¬ 
pute numerically, but there are comparatively few useful analytical solutions (in ID probably 
restricted to the solutions of Riemann’s differential equation [T]). The ease of numerical sim¬ 
ulation is in some respects unfortunate because it leads to the belief that everything is well 
understood about linear wave equations. Finding solutions may be numerically trivial, but 
it is almost impossible to make any general statements on this basis. In this respect trans¬ 
formation optics las] has been a big step forward, providing a simple intuitive recipe for 
relating inhomogeneous materials to analytical solutions (revealing hidden symmetries of 
apparently unsymmetrical structures 0). It is worth keeping in mind that concealment 
devices only became a concrete possibility on the basis of this mathematical insight 

that coordinate deformations are equivalent to inhomogeneous media. 

Aspects of transformation optics have been known for some time: Tamm explored the 
connection between wave propagation in anisotropic media and curved geometries in the 
1920s [7]. The perfectly matched layer [H] is another familiar example, and it was established 
some time ago that a complex transformation of the spatial coordinates is equivalent to an 
inhomogeneous non-reflecting absorbing medium p]. More recently there has been a revived 
interest in the idea of using complex coordinates to understand wave propagation [T0HT2], as 
well as the discovery of a relationship between media with loss and gain and invisibility |T3T 
[15], and an improvement in the practical realisation of materials with tailored dispersion and 
dissipation [12]. Our work is similar in spirit, as we also make use of complex coordinates 
to understand the effect of inhomogeneous material properties on the propagation of waves. 
However, rather than use complex coordinate transformations to derive different sets of 
material properties, we seek to understand the behaviour of the wave in the entire complex 
position plane, applying this as a method for understanding wave propagation in a family 
of inhomogeneous media. 

The key message is that it can be useful to analytically continue a wave equation into 
complex coordinates, and that we can for instance hnd families of reflectionless inhomoge¬ 
neous media through doing this. The particular wave equation we shall use to demonstrate 
this is the equation for the electric held of a monochromatic (frequency u) TE polarized 
wave z ■ E = ip propagating through an isotropic inhomogeneous slab of material 

where propagation occurs in the x-y plane, ko = oj/c, ky determines the angle of incidence, 
and the permeability and permittivity depend on x (see hgure[^. Equation ([^ also holds 
for the TM polarization, with the roles of e and /i interchanged. Although we have written 
([^ as an equation for the electromagnetic held, the same equation occurs in other areas 
of wave physics. For example an acoustic wave obeys ([^, but the bulk modulus and the 
density play the role of the permittivity and permeability. 


ip{x) = 0 ( 1 ) 


I. A PHYSICAL MEANING FOR COMPLEX COORDINATES 

Although ([^ holds for a wave propagating along the real axis x G (— 00 , 00 ), we can 
analytically continue ip into the complex position plane z = Xi + 1 x 2 - The immediate 
thought might be that this is nothing more than mathematical pontihcation. Yet there is a 
well dehned physical meaning for the whole complex position plane, which we now describe. 
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FIG. 1: The solution to the monochromatic wave equation 0 represents a TE polarized wave 
propagating at an angle (determined by ky) through an inhomogeneous medium defined in terms of 
the permittivity and permeability, e(x) and /r(x). In this work we examine the analytic continuation 
of this wave propagation into the complex x plane, and through doing this find profiles e(x) (/i = 1) 
that do not reflect radiation, in certain cases independent of the angle of incidence. 

Writing 0 along a trajectory in the complex plane ^ = z{'^) parametrized by a real 
number 7 the equation becomes 

where z' = dz/d'y. If we interpret 7 as a new position variable, this is equivalent to propa¬ 
gation through a new inhomogeneous anisotropic medium where 

ezzi'y) = e(7)^'(7) 
hxx(7) = h(7)(^'(7))"^ 

hro(7) = h(7)^'(7)- (3) 

In general these material parameters will be complex functions of position, and the medium 
will exhibit some combination of dissipation and gain. The equivalence (0 has an interesting 
consequence: having solved the wave equation as a function of x, analytic continuation 
into the complex position plane is equivalent to solving the wave equation in an inhnite 
number of closely related inhomogeneous media, one for every trajectory ^( 7 ). The material 
parameters will be generally anisotropic, unless ^( 7 ) is parallel to the real axis. Figure 
shows the simplest case of such an analytic continuation, (p{x) —)■ (p(z) = exp(ifcz) with 
two examples to illustrate the equivalence of complex trajectories and materials. Note 
that the equivalence is not restricted to this one dimensional example, and also holds for 
the three dimensional wave equation, as is well known by those working on the theory of 


99(7) = 0 (2) 
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perfectly matched layers [9]. The substitution (p —)■ reduces equation Q to the form 

(p" + (fcoCeff — ^y)'-P — 0) where fcoCefr = — (3/4)(/i'//i)^ + (l/2)/i"//i, arid fj! = d^/dz. 

For the remainder of this work we therefore consider —without loss of generality —the case 
/X = 1 in Q, understanding e as eefr- 



Xi Xi 

FIG. 2: Two examples showing the equivalence between wave propagation along trajectories in 
complex coordinates and in inhomogeneous media. In both panels we plot the wave (p(z) = 
exp(27riz/A) as a function of z = Xi+\X 2 with brightness indicating magnitude and colour indicating 
phase: [0,7r/2, vr, 37r/2] represented as [red,green,cyan,purple]. The top subfigure within each panel 
shows if evaluated along the contour ^(y) shown in blue (in each top subfigure blue is the real 
part, green the imaginary part, and red the magnitude), and the bottom subfigure shows z'{'^), 
which gives the equivalent material parameters via ([^ with e(7) = g{y) = 1. (a) The contour 
-2(7) = 7 + i (1 + tanh(37)) represents an inhomogeneous reflectionless absorber; (b) the closed 
contour 2(7) = 3cos(7) + isin(7) represents a parity-time symmetric periodic medium where the 
wave propagates with a phase that advances then reverses through the unit cell.. 


II. POLES IN THE PERMITTIVITY AND BRANCH CUTS IN THE WAVE 


Having provided a physical motivation for the analytic continuation of (p into the complex 
position plane, we free ourselves to investigate the general behaviour of (p as a function of z. 
In particular we shall show how the phenomenon of reflection manifests itself in the complex 
position plane, and how it can be controlled through the function e{z). 

For the sake of simplicity we restrict this investigation to bounded permittivity prohles 
tending to unity at large \z\, and we choose to represent such an e{z) as an inhnite product 


00 


'(7=n 

i=0 


z- qi 

z-Pi 


(4) 
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where the zeros g* and poles pi are all conhned to a region \z\ < |^max|- The behaviour of the 
wave ip in the complex position plane has some generic features that arise from the positions 
and weights of the poles of Q. 

We have found that the effect of the poles in e{z) on the wave p{z) is to introduce branch 
cuts that run from the poles to inhnity. This can be simply understood within the Born 
approximation mi. which we apply along lines of constant X 2 in the z plane (Green function; 
exp(iA;|a;i — x'i\)/2ik). As xi —)■ ±oo a wave incident from the left can be approximated to 


(p(xi + 1x2) ~ 


ifcxig kx2 dx'[e{x' + 1 x 2 ) — Xi —)■ —cxd 

]^^ikxi^-kx2 dx'[e{x' + 1x2) — 1 ] Xi —)■ cxD 


with k = (kQ — kyY^"^. Given that k > 0, the integration contour in the first of the expressions 
(|^ can be closed in the upper half x' plane, which gives a null result if e{z) is analytic in 
the region Im[z] > X 2 , and a non-zero result if e{z) is not analytic in this region. Note 
that there is a subtlety here: if the permittivity approaches unity as 1 /z or slower then the 
waves won’t be plane waves at inhnity, as assumed in ([^ and in these cases ([^ neglects 
a logarithmic term in the phase. Nevertheless the equation correctly reproduces the hrst 
order rehection and transmission coefficients (see appendix [A]) . 

In the simplest case of a single pole €{z) = {z — qo)/{z — po) the hrst of (|^ evaluates to 


(p(xi - 00 ) ~ e’^^ + e-’^" 


0 


Im[po] > Im[ 2 ;] 
Im[po] < Im[ 2 ;] 


( 6 ) 


Meanwhile on the far right of the prohle, xi —)■ cxd 


(p(xi —)■ 00 ) 




ikz 



Im[po] > Im[ 2 ;] 
Im[po] < Im[ 2 ;] 


(7) 


Applying the boundary condition (p(xi —)■ 00 ) ~ exp (ikz), it is clear that (p(z) has a branch 
cut that runs horizontally left from the pole at Zq. If the wave had been incident from the 
right, then the branch cut would run to the right. For the remainder of this work we restrict 
ourselves to the case of incidence from the left, with the understanding that similar results 
hold for the other direction. Figure shows an illustration of this phenomenon, where we 
have numerically integrated the wave equation and conhrmed the presence of the branch 
cut for this simple case. In appendix we give an analytic example, also confirming the 
existence of the branch cut and the validity of ([ 6 ]-[^. 

In general a branch cut in (p(z) runs from each of the poles of e(z) to inhnity, and we 
shall now argue that the presence of the cuts is due to the phenomenon of rehection. Firstly 
we note that in complex analysis, Liouville’s theorem [18] shows that every bounded entire 
function is constant. This means that an inhomogeneous permittivity e(z) that tends to 
a constant value at large \z\ cannot be an entire function and must contain poles. Now 
consider the illustration given in hgure|^ Suppose we move out to the semi-circle |z| —>■ cx) 
in the upper half plane, and consider wave propagation along the contour ^( 7 ) = i?exp (iy). 
Using ([^ we hnd that this is equivalent to propagation through a medium where €zz = 
[1 + 0(l/i?)]ii?exp (iy), pyy = iRexp (iy) & = —iexp (—i'y)/R. This is a material with 

a negative index and a large degree of gain in the region y G [7r/2, tt] and a negative index 
and a large degree of dissipation in the region y G [ 0 , 7 r/ 2 ]. Starting from the boundary 
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FIG. 3 : Numerical solution of Q for the function ip(z) in the profile €(z) = 1 — ‘laji^z — zq), fJ. = 1 , 
for the case ky = 0 , a = 1 , ko = 1 , and zq = i. The numerical solution was obtained through 
imposing the boundary condition (p{z) = 1 and (p'{z) = iy^ e{z)ko at some position on the far right 
of the profile where e'{z) ~ 0, and then numerically integrating Q along a series of lines parallel 
to the real axis. The red dot at zq indicates the pole in the permittivity profile, and the red dashed 
line indicates a branch cut. The two lower panels (i) and (ii) show the functional form of ^{z) 
evaluated along the white dashed lines shown in the upper panel (colour conventions as in figure]^ 
and the inset panels show the spatial variation of e along these lines. The wave in panel (i) is 
evaluated in a region just above the branch cut where the medium is dissipative, and despite the 
rapid variation of e(z) there is no reflection, evident in the lack of oscillations in the magnitude of 
the wave where x < 0. Panel (ii) shows the wave evaluated below the branch cut in a region where 
the medium is amplifying: the oscillations in the absolute value of the field on the left show that 
here there is significant reflection. 

condition that the wave is right-going at 7 = 0, and integrating throngh the dissipative 
region, the wave is exponentially diminished at 7 = 7 i/ 2 . Meanwhile, starting from the 
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FIG. 4 : Imposing the boundary condition that the wave is right-going on the far right of the profile 
= (/?_!_, and integrating the wave equation Q around a large semicircle in the upper half plane 
^(7) = R exp (17), we find that the wave becomes exponentially small. Meanwhile, were we to follow 
the same procedure on the left and follow a sum of left and right-going waves ip = + Bip- 

along 2(7) in the clockwise sense then we would conversely find that the solution was exponentially 
large in the upper half plane. The presence of branch cuts (red dashed lines emerging from the 
poles of e{z)) in (p avoids this contradiction. The wave is analytic in the region of the complex 
position plane X2 > l^maxl) where it also tends to zero as 2:2 — >• 00. All the inhomogeneous media 
captured in this region are reflectionless. 


boundary condition of a sum of right and left going waves at xi = —00, and integrating 
through the gain medium, the wave is exponentially amplihed at 7 = 7r/2, unless 5 = 0 . 
This leads to a contradiction unl ess there is a branch cut in ip that forces the reflected wave 
to disappear past a certain value of 7. This is the jump captured in expression ([^. Note 
that the branch cut can be moved to run from the pole to inhnity in any way we choose, 
although it is non-physical to consider wave propagation along a contour that passes across 
the cut. Moreover this shows that above all of the cuts the reflection disappears, even though 
the medium may be rapidly varying in this region of the complex position plane. 

We can use this understanding to hnd a family of inhomogeneous media that reflect no 
radiation for any angle of incidence—a hnding already reported in IZB, but derived in a 
different way. If we consider the region of complex position above all of the poles in e(z), 
^2 > kmaxi (see hgure and the branch cuts in p run parallel to the Xi axis, then in that 
region p{z) is an analytic function that tends to zero as 0:2 —)■ C)0. In analogy to what is 
typically done in the frequency domain | 22 ] , such a function can be represented as a Fourier 
integral over positive wave-numbers 

dk 

Jo 27 r 

The lack of any negative wave-numbers means that the reflection is zero for incidence from 
the left of all the prohles captured in this region, and this is independent of the value 
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of ky (because ky doesn’t change the position of the poles in ([^). This means that all 
permittivity profiles that are analytic at complex positions above the propagation axis do not 
reflect radiation incident from the left. Figures]^ and|^i conhrm this fact for the simple case 
of a single pole in the permittivity prohle. Figure illustrates wave propagation above the 
pole at Zq where the reflection is zero, and £gure|^i shows propagation below the pole, where 
the permittivity varies less rapidly but strong reflection is evident. Appendixcontains an 
analytic demonstration of this phenomenon in a simple exactly solvable case. 


III. NO BRANCH CUTS, NO REFLECTION 


Given that the branch cuts in tp are intimately connected to the phenomenon of reflection, 
we now investigate under what circumstances they disappear and thereby determine a large 
set of inhomogeneous media from which the reflection is zero (remembering that the complex 
position plane describes wave propagation in a family of inhomogeneous media). To avoid 
confusion we must emphasise that while the disappearance of the branch cuts is sufficient 
to remove reflections, it is not necessary: we know that if there is reflection then there will 
be a branch cut, but this doesn’t imply that if there is a branch cut then there is reflection. 
We assume an ansatz for the wave (p that is free of branch cuts and therefore by dehnition 
free from reflection, 

^{z) = F{zy^^ ( 8 ) 


The function F[z) is assumed to be without branch cuts, with zeros at isolated points r* 
and poles at points Sj. Assuming that F —>■ 1 as | 2 ;| —)■ cxd, consistent with the form of e{z) 
given in (|^, we write 


M 


m = n 


i=0 


(z — ri)"^^ 

{Z - Si)'^i ' 


(9) 


where mi and ni are integers such that the numerator and denominator of ([^ are polynomials 
of the same degree, and M is an integer which can be formally taken to inhnity. Inserting 
(|^ in ([^ with /i = 1, we hnd that for such a form of ip the permittivity must equal 


1 - P 
1 


F(z) 

dG{z] 

dz 


+ 2ik 


F'i 


F(z) 

+ lG{z)Y + 2ikG{z) 


where F' = GF, with 


= E 


mi 


ni 


ri 


Si 


( 10 ) 


( 11 ) 


Therefore, after constructing a function G{z) as the series of simple poles 0 with 
T.^1 = '^ni, the permittivity prohle ( [Io| ) will be such that p is free from branch cuts 
and the rehection vanishes for the whole complex position plane. However, in general this 
is dependent on the value of ky because (10) explicitly depends on k. We now show how 


this angle dependence can be eliminated to give zero rehection for all angles of incidence. 
Separating out the expression into a sum of hrst and second order poles, the expression for 















the permittivity becomes 


e(z) 



mi{mi - 1) ni{ni + 1) a; k 

{z-ri)‘^ (z-SiY z-ri z - Si 


( 12 ) 


where 


ai = 2mi 

bi = 2ni 


qr-, _ 7 ^ ly j _ 


p^i 


n-Tp ^ ri- Sp 


-ifc+y^^ + y—^ 


Si Sp ^ Tp Si 


(13) 


This tells us how to relate the residues of the poles in the permittivity prohle to eliminate the 
branch cuts discussed in the previous section, which is numerically demonstrated in hgure 
As already discussed, in general (12) is a function of ky, because it depends on k through a; 


and bi. This means that a different inhomogeneous medium is required to suppress reflection 
for each angle of incidence. To eliminate reflection for all angles of incidence requires this ky 


dependence to disappear from (13), which for example would require mi = 1 and a choice 


of ri such that ai = 0, and that bi is independent of k. 

The simplest example of (12) is to take a single pole and zero in (|^ with mi = ni = 1. 
This leads to the following family of reflectionless prohles 


^0 


(z - Si)2 


+ 


ifc(ri — Si) — 1 

- Si) 


z — ri z — Si 


(14) 


which are complex functions of position along the lines of constant X 2 (see hgure In the 


)-2 


particular case where ri = si—i/k we obtain the permittivity prohle e{z) = l — 2kQ^z—si) 
which is rehectionless for all ky. Similarly, taking ni = n, nj>i = 0, mi<n = 1 and = 0 
we hnd the permittivity prohle 


1 - p 


n(n + 1) 
(^-si)2 


ai 


'y. — 


ri 


+ 


bi 


z — Si 


where ai = 2[ik + ~ ''"p) ~ '^1 ~ -Si)] and bi = — The n quantities ri can 

be chosen so that a; = 0 for all /, and then bi is automatically zero. One therefore hnds 
that e{z) = 1 — n(n + l)/C(C^(z — Si)“^ is rehectionless for all angles of incidence when n is 
an integer. These omni-directional rehectionless prohles were investigated some time ago 
by Berry and Howls [12] in their considerations of the WKB approximation applied to the 
Poschl-Teller potential [20], and they are also a special case of the functions discussed in 
appendix 

The aforementioned rehectionless complex prohle is actually a constitutive element of 
the Poschl-Teller prohle 1271 , which is also a special case of ([l2|-[l^. Taking ruj = 1, 


Sin = (I + l/2)i7ra and = n (all other s, and rii are zero) we have 


c(z) — 1 + 


n(n 


1 ) 


kQa'^cosh'^ {z / a) k^ 


OO 


l=—oo 


n—1 

E 

_p=0 




+ 


Ml 


Tnl+p Z-i{l+ 1)™ 
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FIG. 5: Complex plane representation of wave propagation in the permittivity profile €(z) = 

where o = [i/c(ri — si) — l]/(ri — si), with n = 5/2 — i and 
Si = 1 + i. This permittivity reduces to (14) when A = B = 1. Panels (a-c) show that as A and 


1 2 

^0 L 


B 

{z-siY 


+ CL 




Z — Sl 


B are both brought towards unity the branch cuts emanating from the poles in the permittivity 
at ri and si disappear. The insets of each of these panels show the permittivity evaluated along 
the real line. Panels (d-f) show the wave propagation evaluated along the real axis in each of these 
cases, and they demonstrate that once the branch cuts have disappeared (strictly speaking only 
the branch cut from si has to disappear), the reflection from the profile vanishes. Note that in this 
particular case only a very slight change to the form of the permittivity profile between (b) and 
(c) leads to a large change in the reflection in panels (e) and (f). 

Assuming rin+p = {I + l/2)i7ra + the coefficients ai and bi become 


Otln+p 2 


n—1 


ik -coth 

n 


ti=0 




bln = 2n 


n—1 


-ik H— coth 


u=0 


a 


^coth 
a \ a 


n—1 

p=0 


(15) 
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In the same way as previously discussed, the a* can be chosen to make the a* zero. Having 
done this the bi are also automatically zero. Thus from the assumption of the absence of 
branch cuts in the complex position plane we reproduce the result that the Poschl-Teller 
prohle e{x) = 1 + n{n + l)(A;oa)~^ cosh^(x/a) is reflectionless for all ky. Indeed, the above 
analysis does not rely on any assumptions on the value of a, nor on the choice of origin of 
the complex position plane. Therefore if we choose a = |a| exp(i6') as any complex number, 
and replace z \yY z — where zq = l^^ol exp(i0), then the Poschl-Teller potential remains 
reflectionless for all angles of incidence 


6(Z) 


1 + 


e + 1) 

fcolaP cosh^[( 2 ; — \zo\e^'f')e~'^^/\a\] 


(16) 


a fact that cannot be established on the basis of the inverse scattering method mm, 
because it assumes real values for the permittivity. 



FIG. 6: The Poschl-Teller potential (16) is reflectionless for all angles when n is an integer. The 
profile is made up of a line of poles parallel to the imaginary position axis, (a) shows wave 
propagation in this potential when re = 1, demonstrating the lack of any branch cuts in the wave 
when the profile is reflectionless. Meanwhile (b) shows wave propagation when re is non-integer, 
showing that branch cuts emerge from each of the poles. 


IV. RELATION TO THE WKB METHOD 

The phase integral (WKB) method [23] already uses an analytic continuation of the wave 
equation into complex coordinates. At first sight the findings of this paper seem at odds 
with the known results of this method, which emphasize that the zeros of e{z) rather than 
the poles are most important for determining the reflection. In this section we show how to 
interpret our findings in terms of the WKB solutions to the wave equation. 

For now, consider the case where ky = 0. The WKB approximations to the solutions of 
the wave equation ([^ may be found in standard texts on phase integral methods, e.g. m 


11 














and are 


1 


(17) 




e(z 


:.±ifeo y/e(z')dz' 


These expressions are valid approximations when fco 3> i.e. when the permittivity 

varies signihcantly on a scale mnch larger than the wavelength (in particnlar, these expres¬ 
sions are an excellent approximation for large \z\, where e varies slowly and is close to nnity). 
The reference point Zr is where the phase of the wave is zero, and we choose this to be a 
complex position at which e(zr) = 0. For the sake of consistency with the literatnre |23] we 
will denote the two expressions in ( [ItI ) by = {zr, z) and = {z, z^)- 

In the WKB method, the phenomenon of reflection is associated with the breakdown of 
the approximation at the zeros of e{z), and the complex positions of these zeros are often 
nsed to calcnlate reflection coefficients [25]. In order to find the correct WKB approximation 
to a particnlar exact solntion to the wave eqnation, we mnst nse a patchwork of different 
linear combinations of the two WKB approximations ip = A{zr, z) +B{z, Zr) thronghont the 
complex position plane. Ultimately it is the changes in these coefficients A and B as we 
move throngh the complex position plane that determines the amonnt of reflection from the 
medinm. The change in these coefficients is known as the Stokes phenomenon^ and occnrs 
across what are known as Stokes lines, defined as the cnrves in the complex plane satisfying 


Re 



= 0 


(18) 


complementary to these lines are the anti-Stokes lines, defined as the cnrves satisfying 


Im 



= 0 . 


(19) 


On the anti-Stokes lines both {zr, z) and {z, Zr) have the same magnitnde, becanse the phase 
is pnrely real. These cnrves divide the complex plane into regions where one of the two WKB 
solntions has a larger amplitnde than the other. Typically the larger of the two is called the 
dominant solntion, and the smaller is called the sub-dominant solntion. The Stokes lines are 
the cnrves along which the amplitnde of the dominant WKB solntion is maximal compared 
to the snb-dominant solntion. In general as we cross a Stokes line the error inherent in 
the WKB expression for the dominant wave will exceed the magnitnde of the snbdominant 
wave. In this region the coefficient of the snbdominant wave is nndetermined, and in general 
it mnst be changed after one has passed throngh the Stokes line in order that the patchwork 
of WKB solntions represent a good approximation to the exact solntion. In WKB theory, 
reflection is associated with this change in the coefficient of the snbdominant wave across 
a Stokes line, as illustrated in figure This is the Stokes phenomenon. As \z\ —>■ oo, 
e{z) —?• 1 and we can see from the definitions (18-19) that the Stokes and anti-Stokes lines 
will asymptote to lines parallel to the imaginary and real position axis, respectively. 

To compare with the findings of the previous sections, we now consider a permittivity 
prohle containing poles only in the lower half position plane, which we have suggested ought 
to be reflectionless for all angles of incidence. In figure]^ we consider a purely right travelling 
wave on the far right (xi —)■ cx)) of a generic prohle containing poles in the lower half position 
plane. Moving along a large semi-circle in the upper half plane (where a right travelling wave 
becomes snbdominant before encountering any Stokes lines), we hnd that the conhguration 
of Stokes and anti-Stokes lines as |; 2 | —)■ cx) guarantees zero rehection |28|. Notice that the 
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FIG. 7: The Stokes and anti-Stokes lines for the two permittivity profiles (a) e{z) = 1 — ^ 
and (b) e{z) = 1 — ^ 2 (subscripts ‘s’ and ‘d’ indicate the relative dominancy of the two WKB 

(-S+l) 

solutions, and insets (i) and (ii) show the permittivity profiles along the real line, with the same 
colour conventions as in the previous figures). Both panels illustrate change in a general WKB 
approximation taken from the right anti-clockwise about Zr across the Stokes line. Both profiles 
show the same behaviour illustrated in figure The small amplitude right propagating wave 
is unchanged across this curve, whereas the large amplitude left propagating wave undergoes a 
change which depends on a constant T, called the Stokes constant associated with the line. For 
completeness we have included the branch cuts that occur in the WKB approximations ( [l7| ), which 
are purely due to the approximate nature of these expressions [2Tj • 


branch cuts emerging from the poles that we discussed previously must be incorporated into 
the WKB theory if we are to avoid the conclusion that the medium is also reflectionless for 
waves incident from the right (see |21j). 

The WKB method is therefore consistent with the previous sections, where we found that 
the region above the poles in e{z) corresponds to a family of reflectionless inhomogeneous 
media for radiation incident from the left. Moreover, via this method we can also gain 
further information about the transmission through these prohles. When considering a 
complex valued permittivity the position of the reference point Zr in ( [ItI ) is not arbitrary. 
This is because both amplitude and phase change with Zr- We take the reference point such 
that a right-going wave has unit amplitude on the far left of the prohle Xi —)■ —00 (shown 
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— Anti-Stokes 
— Stokes 

X\— Branch cut in exact solution 
WKB branch cut 


FIG. 8: A schematic diagram of the complex position plane for a profile analytic in the upper half 
plane, containing poles at pi and zeros at qj. Outside the dashed black circle, the permittivity is 
approximately one, so the Stokes and anti-Stokes lines are approximately horizontal and vertical, 
respectively. When the radius of the semi-circle is sufficiently large, we will cross a number of anti- 
Stokes lines, followed by a number of Stokes lines and then more anti-Stokes lines before reaching 
the negative real axis. Just before encountering the Stokes lines, {zr.,z) will be a subdominant 
(low amplitude wave). Because of this subdominancy, crossing the Stokes lines will not change the 
wave and we arrive at the negative real position axis with the solution remaining solely as a right 
travelling wave. Hence we get zero reflection, in agreement with the argument of section |IT| 


as point b in figure]^. The amplitude of the wave on the far right of the prohle (point a of 
hgurej^ then gives us the transmission coefficient, 

t= lim _ .20) 

a^oo ^ ' 

b^—oo 


Note that the argument of this quantity will not converge, so we cannot strictly give a 
meaning to the phase of the transmission coefficient through such infinitely extended profiles. 
However, the time average of the transmitted power is meaningful, which is 


|2 _ -2fcoIm(/p 


( 21 ) 


where the integration contour has been deformed into a clockwise semi-circle in the upper 
half position plane. Note that for a permittivity profile composed of a collection of poles in 
the lower half plane, e{z) = 1 -|- Q!i/(^ — ^i) «„/{z — Zn) + f3i/{z — Zn+i)"^ + ■ ■ ■ + 


/3nl{z — Zn+m) + • • •) (21) depends only on the sum of the residues of the simple poles 


|f|2 = 


Therefore, for a permittivity profile that is analytic above the line of propagation, in addition 
to zero reflection for incidence from the left, the transmission is unity when the sum of the 
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residues of the simple poles of e{z) is zero. The reflection coefficients are less trivial to hnd 
because they are dependent on the values of the Stokes constants, which in general we do 
not know. 


V. SUMMARY AND CONCLUSIONS 

In this work we have explored the utility of analytically continuing wave equations into 
the complex position plane. In a similar spirit to the formalism of transformation optics, 
we found that in electromagnetism there is a simple interpretation for the entire complex 
position plane, where each curve ^( 7 ) can be understood as a different inhomogeneous 
anisotropic medium. We then investigated the properties of solutions to the one dimensional 
Helmholtz equation ip" + [k^e^z) — ky](p = 0 in the whole complex position plane, using a 
combination of exact and approximate analytical techniques, including the WKB method. 

Considering bounded permittivity prohles constructed as a sum of poles of varying de¬ 
grees, we found that in general there are branch cuts in the wave p>{z) that emanate from 
the poles of e{z). These branch cuts are connected to the reflection of the wave. In the 
region of the complex position plane that lies above the poles in the permittivity there is 
no reflection, whereas in general below there is, and the branch cuts account for this jump. 
We found that this knowledge can be used to construct reflectionless permittivity prohles in 
two distinct ways. We can either construct a prohle where all the poles occur in the region 
of complex position below the axis of propagation (which reproduces the recent hndings 
of [ 21 ]), or we can demand that the branch cuts disappear, which we have shown reproduces 
a generalisation of the Poschl-Teller potential, equivalent to a permittivity prohle in optics. 
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Appendix A: An analytic solution 

In this appendix we conhrm—using a particular family of exact solutions to ([^— that 
branch cuts in the wave emerge from poles in the permittivity. The particular equation of 
interest is a type of conhuent hypergeometric equation: Whittaker’s diherential equation [1] , 
which takes the form 




and has two standard solutions, H±K,/i(±C)- Changing variables ( = 2ikQ{z + izq) we hnd a 
diherential equation in the form of (j^ with /r = 1 , = 0 and a permittivity 



which has a pole at the complex position —i^o- Close to —izo the solutions take a particularly 
simple form where the branch cut can be seen immediately. When /i 7 ^ 1/2, the double pole 
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dominates and the solution approximates to 


(p ~ (z — 

and when p = 1/2, the simple pole dominates and the solution is 

^ + {z- iZo) In (^ 0 ( 2 : - i^o)) 




{z - i^o) - Ik/co [z - izoY 


(A2) 


(A3) 


Notice that in general both (A2) and (A3) have branch cuts emerging from —izQ. These are 
the branch cuts discussed in the main text. Now we examine the opposite case, and hnd the 
behaviour of these branch cuts far from the pole. 

The exact solutions 1 T±k,^(±(/)(z) are cut along the line arg(z) = ivr, and the general 
formula to continue the wave across this cut onto the next Riemann sheet is given in terms 
of the Whittaker functions evaluated on the hrst sheet by [1] 

(-irw^jze^-n = + 8m((2m - 2),z7r) 


sin(2p,7r) 


sin(2m/i7r)27rie'’ 


(2p7r)r (^ + p - k) T (I - /i - k) 


sm 




(A4) 


Now consider the form of these functions along the real z axis at ±cx). To hnd this we use 
the asymptotic forms of the Whittaker functions on the hrst sheet 


W^^^{\z\ -)■ 00 ) 




(AS) 


We take the particular case of /i = | and k = For the solution kF_i^ i{—2ikQ{z + izQ)), 
with zq > 0 the argument remains on the hrst sheet as we trace the solution from z = +00 
to z = — 00 . Therefore asymptotically the wave is 


W i. 


1 \ 

2 ’2 


-2iko(z + izo)) 


^iko (z+i^o) g- [log( 2 fco |z+izo |) - ] 

, (z+izQ) p,- [log( 2 fco 1 2+izo I) - ^] 


^ —)■ +CXD 

^ —>■ —00 


(A6) 


e 2 g‘"-uv-^'-u;g 

which is right-going on both sides and exhibits no rehection. Meanwhile if < 0, we pass 
through the branch cut as we follow the solution back to z — 00 . Applying (A4) we can 
move the cut parallel to the real z axis and we hnd the asymptotic forms 

i(-2iko(z + izo)) ~ ^ ^ +00 (A7) 


and 


kF_i^ 1 (-2ifco(^ + Ro)) ~ ] 


27ri 


r(i^)r(i + k|T) 


g- ifco (2+120) g [log( 2 fco 1 2+izo I) + T 


^ z —)■ —00 (A8) 


2 V-" ' 2 

As stated in the main text, the presence of the cut is clearly connected to the phenomenon 
of rehection. As we move from above to below the pole at —izo, the solution to (Al) goes 
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from ( |A 6 [ ) to ( |A7 -A8), and on the left we now have a reflected wave. The above analysis 
enables to identify the reflection and transmission coefficients of this extended profile, which 
are (above the cut) 


T>=e- 
= 0 


/cqttA 

2 


(A9) 


and (below the cut) 


= e 2 


R< = 


27re 


/cqttA 


r(i^)r(i + had) 


fcnTrA 

= 2 e 2 


smh(-) 


(AlO) 


where to obtain the second of ( |A 10 [ ) we applied the properties of the T function given 
in [T] (eq. 5.4.3). To first order in A these reflection and transmission coefficients are 
T>/< ~ 1 ^ ^ and i? 
given in 


~ i -|- auu rt< ~ nkoA, in agreement with the results of the Born approximation 


In cases where k = 0 and /i 7 ^ 0 in (Al), the permittivity approaches unity asl /and the 
reflection coefficient i?< has a denominator containing r(^ ± /i). When fi equals an integer 
plus one half we are at a pole of the T function, and the reflection coefficient vanishes. In 
this case the branch cut also vanishes, and we reproduce the reflectionless behaviour found 
in section m 
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